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Abstract. We combine samples of spiral galaxies and starburst systems observed with ISOCAM on board ISO to 
investigate the reliability of mid-infrared dust emission as a quantitative tracer of star formation activity. The total 
sample covers very diverse galactic environments and probes a much wider dynamic range in star formation rate 
density than previous similar studies. We find that both the monochromatic 15 /im continuum and the 5 — 8.5 /im 
emission constitute excellent indicators of the star formation rate as quantified by the Lyman continuum luminosity 
Lhyc, within specified validity limits which are different for the two tracers. Normalized to projected surface area, 
the 15 pm continuum luminosity £i 5Mm , ct is directly proportional to El Y c over several orders of magnitude. Two 
regimes are distinguished from the relative offsets in the observed relationship: the proportionality factor increases 
by a factor of fa 5 between quiescent disks in spiral galaxies, and moderate to extreme star-forming environments 
in circumnuclear regions of spirals and in starburst systems. The transition occurs near El Y c ~ 10 2 Lq pc -2 and 
is interpreted as due to very small dust grains starting to dominate the emission at 15 /jm over aromatic species 
above this threshold. The 5 — 8.5 )im luminosity per unit projected area is also directly proportional to the Lyman 
continuum luminosity, with a single conversion factor from the most quiescent objects included in the sample up 
to Sl Y c ~ 10 4 Lq pc -2 , where the relationship then flattens. The turnover is attributed to depletion of aromatic 
band carriers in the harsher conditions prevailing in extreme starburst environments. The observed relationships 
provide empirical calibrations useful for estimating star formation rates from mid-infrared observations, much 
less affected by extinction than optical and near-infrared tracers in deeply embedded Hn regions and obscured 
starbursts, as well as for theoretical predictions from evolutionary synthesis models. 
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1. Introduction 

Star formation is a fundamental process of galaxy forma- 
tion and evolution. Estimates of the star formation rate 
(SFR) in galaxies at all redshifts are key indicators of the 
efficiency and mechanical feedback effects of star forma- 
tion activity, of the chemical evolution of the interstellar 
and intergalactic medium, and, ultimately, of the cosmic 
star formation history. 

Commonly used probes of the SFR include photo- 
spheric emission from hot stars in the ultraviolet, nebular 
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H and He recombination lines as well as fine-structure lines 
arising in Hn regions from optical to radio wavelengths, 
and the total infrared luminosity (A = 8 — 1000 /xm), 
the bulk of which is due to heated dust reprocessing 
the interstellar radiation field (see, e.g., the review by 
lKennicuttlll998|) . However, ultraviolet, optical, and even 
near-infrared diagnostics are subject to potentially large 
uncertainties because of extinction in deeply embedded 
young star-forming sites and in nuclear regions of galaxies. 
Nebular lines may be difficult to measure when intrinsi- 
cally weak or superposed over a strong continuum. While 
dust emission suffers very little from extinction effects and 
is usually strong in star-forming environments, the total 
infrared luminosity is difficult to evaluate because it is 
generally derived from observations in a few wavelength 
intervals which do not constrain the spectral energy dis- 
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tribution accurately. Moreover, a cirrus-like component 
unrelated to star-forming regions can contrib ute substan- 
tiall y to the far-infrared output of galaxies l|Heloulll98fil 
ISauvaee fe Thuanlll992j) . 

Mid-infrared emission (MIR, A = 5 — 20 /im) provides 
an alternative probe of star formation activity. The "clas- 
sical" spectrum of star-forming sources exhibits broad 
emission features often referred to as "unidentified in- 
frared bands" (UIBs) and of which the most prominent 
dominate the 6 — 13 /im range, and a continuum rising to- 
wa rds long wavel e ngths at A ^ 11 pm (see the reviews 
bv iGeballd Il997t iTokunagal Il997t Icesarskv fe Sauvaeel 
Il999t iGenzel fc Cesarskvl 1200011 . Various carbonaceous 
materials have been proposed to carry the UIBs, in- 
cluding the popular polycyclic a romatic hydrocarbons 
CPAH s: e.g. |Leger fc Pugetl Il984h that we adopt here- 
after. IPeeters et alJ l)2002f) have analysed their shape 
and relative amplitude variations in different classes of 
Galactic objects. The continuum emission is generally at- 
tribu ted to very small dust grains (VSGs; e.g. lDesert et, alJ 
1990) about which little is known. Superposed on these 
PAH and VSG components, H recombination lines and 
fine-structure lines of various metals originating in H n 
and photodissocia tion regions are observed as well (e.g. 
ISturm et al J 1200(1) . These lines may however not always 
be measurable because of their weakness or of insufficient 
spectral resolution. 

Numerous past studies have established that PAH and 
A <; 1 1 /im continuum emission trace well star- forming 
regions but their usefulness as quantitative diagnostics 
is still debated. Complications arise from the different 
nature of the emitting particles and by their being out 
of thermal equilibrium under most radiation field condi- 
tions, unde rgoing; large temperature f l uctuations of several 
100 K fe.g.lGreenherg fc HoT^IT^t Inline fc AndersonI 
Il98,4 IPuget fc Legerlll989l) . In addition, although both 
species are predominantly heated by energetic radiation, 
PAHs can also be excited by softer optical and near- 
ultraviolet photons as indicated by their detection in the 
diffuse interstellar medium, in regions of insufficient far- 
ultr aviolet energy density to account for their heating, 
fe.g.lSellgren et al.ll99"o]lMattila et alll996tl"Uchida et all 
ll998ll2000HLi fc Drainel2002j) . Furthermore empirical ev- 
idence indicates that the A ^ 11 /im emission is produced 
by a mixture of dust particles akin to PAHs (or at least 
whose flux variations follow well those of PAHs) a nd of 
VSGs (e.g. iHonv et allEoOll IRoussel et a,lJl2001bl). The 
first component is best seen in quiescent environments 
such as disks of spiral galaxies while the second becomes 
prominent in active star formation sites. It is not yet clear 
how their combined emission varies over a large dynamic 
range in star formation intensity. 

On the other hand, spatially resolved studies of 
Galactic and Magellanic Clouds H n regions have revealed 
that both the PAH features and the VSG continuum 
are produced in the vicinity of massive stars, the for- 
mer arising in photodissociation regions (PDRs) at the 
interface between ionized and molecular gas and the lat- 



ter peaking closer to the ionizing stars (e.g.lGeballd 
Tokunagalll997[ IVerstraete et alJll996l ICrete etall 



Contursi et all 



1997 



1999; 



2000). MIR imaging of external galaxies 



has shown that bright emission from both components 
is closely associated with active star-forming sites on 
large scales as well (e.g.lMirabe l et al.lll998HMattila et all 
Il999t IRoussel et alJEoOlcUForster Schreiber et alJl2003j) . 
A strong coupling with the SFR may thus exist and 
has been demonstrat ed for disks of spiral galaxies by 
lR.oussel et alJ l)200ld) . Specifically, these authors found 
that the broadband 5 — 8.5 /mi and 12 — 18 /im fluxes vary 
linearly with the Hct line flux in the disks of 44 spirals. 
lF5rster Schreiber et alJ $2003) also found a direct propor- 
tionality between the monochromatic 15 /mi continuum 
(AA = 0.4 /im) and the [Ar n] 6.99 /im line emission in the 
nearby starbursts M82, NGC253, and NGC 1808, down 

to spatial scales of ~ 100 pc. 

In this paper, we pursue the work of IRoussel et all 
l)2001d) and iForster Schreiber et all l(2003h by combining 
samples of spiral and starburst galaxies obs erved with the 
ISOCAM instrument l|Cesarskv et allll996h on board th e 
Infrared Satellite Observatory (ISO: iKessler et alJll996ft . 
The merged sample covers diverse environments ranging 
from quiescent galactic disks to infrared-luminous merg- 
ing systems. This allows us to extend the investigation to 
higher activity levels and to test whether previous results 
restricted to specific environments can be generalized into 
more universal relationships. We derive the dependence 
of the PAH-dominated 5 — 8.5 /im emission and the VSG- 
probing monochromatic 15 /mi continuum on the produc- 
tion rate of Lyman continuum photons Qhy C quantifying 
the SFR. The resulting empirical calibrations provide use- 
ful tools in MIR studies of star-forming galaxies as well as 
constraints for models predicting the dust emission of such 
systems. 

The paper is organized as follows. Section presents 
the galaxy sample. Section|21describes the MIR indicators 
and the SFR estimates obtained from more classical di- 
agnostics. Section 0] discusses the derived calibrations and 
Sect. summarizes the results. 



2. Galaxy sample 

We drew our sample from separate studies published 
by us and from the ISO archive. All sources were ob- 
served with ISOCAM either with the broadband filters 
LW2 centered at 7 /im (5 — 8.5 /im) and LW3 centered 
at 15 /im (12 — 18 /im) or with the continuously variable 
filter (CVF) covering the 5 — 17 /im range at a resolu- 
tion of R = A/AA ~ 40. Ten galaxies were observed 
in both broadband photometric mode and spectrophoto- 
metric mode. We used this subset to assess the photo- 
metric consistency and to derive conversion factors be- 
tween measurements obtained through the ISOCAM fil- 
ters and from the CVF spectra (Sect. CUD • Details of the 
observations are given in the relevant references (Table^). 
To our knowledge, the ISOCAM data of IC 342 have not 
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Table 1. Sample of galaxies and ISOCAM observations: circumnuclear regions of spiral galaxies and starbursts 



Source 



Distance a Morph. type 6 Nuclear type c 



Observations e Project 



T 



Reference 3 





(Mpc) 






(Lq) 








Circumnuclear regions of spiral galaxies 












NGC986 


23.2 


SBab 


Hn 


4.6 x 10 10 


LW 


Sf glx 


1 3 


NGC 1097 


14.5 


SBb 


Liner /Sy 


3.8 x 10 10 


LW + CVF 


CAMbarre 


2 


NGC 1365 


16.9 


SBb 


Sy 


8.7 x 10 10 


LW + CVF 


CAMspir 


2 


NGC 4102 


17.0 


SABb 


Hn 


4.2 x 10 10 


LW 


SLglx 


1, 3 


NGC 4293 


17.0 


SBO/a 


Liner 


4.7 x 10 9 


LW 


Virgo 


4, 2 


NGC 4691 


22.5 


SBO/a 


Hn 


2.4 x 10 10 


LW 


CAMbarre 


2 


NGC 5194 (M51) 


7.7 


SAbc, int. 


Sy 


2.4 x 10 10 


LW + CVF 


CAMspir 


5, 2 


NGC 5236 (M83) 


4.7 


SABc 


Hn 


1.9 x 10 10 


LW + CVF 


CAMspir 


2 


NGC 6946 


5.5 


SABcd 


Hn 


1.4 x 10 10 


LW + CVF 


SLglx, Zzcvfcam 


6, 2 


NGC 7552 


19.5 


SBab 


Hn 


8.6 x 10 10 


LW 


CAMbarre 


2 


NGC7771f 


57.2 


SBa, int. 


Hn 


2.1 x 10 11 


LW 


SLglx 


1, 3 


Starburst galaxies 
















NGC 253 


2.5 


SABc 


Hn 


1.8 x 10 10 


CVF 


CAMACTIV 


7 


NGC 520 


27.8 


Pec, merger 


Hn 


6.5 x 10 10 


LW + CVF 


CAMACTIV 


8 


NGC 1808 


10.8 


RSABa, int. 


Hn 


3.8 x 10 10 


LW + CVF 


CAMACTIV 


7 


NGC 3034 (M82) 


3.3 


10, int. 


Hn 


4.8 x 10 10 


CVF 


CAMACTIV 


7 


IC 342 


3.3 


SAB(rs)cd 


Hn 


3.3 x 10 9 


CVF 


IMSP.SBG 




LIRGs/ULIRGs 
NGC 3256 


37.4 


Pec, merger 


Hn 


4.0 x 10 11 


CVF 


CAMACTIV 


8, 10 


NGC 6240 


97.2 


10: pec, merger 


Liner 


6.0 x 10 11 


LW + CVF 


CAMACTIV 


8, 10 


IRAS 23128-5919 


180 


merger 


Hn 


9.4 x 10 11 


LW + CVF 


CAMACTIV 


9 


Arp 220 


72.5 


S? (Pec), merger 


Liner 


1.3 x 10 12 


LW + CVF 


CAMACTIV 


8, 10 



a Dist ances are from the NGC ca talogue l|Tullvlll988l) with the following exceptions: NGC 253: fPavidge fc Pritched dl990T) : 
M 82: lFreedman fc Madorel ll 19881) : IC342: lsaha et all (I2002T) : NGC 7771, NGC 624 and Arp 220: computed fro m the Hi red- 
shifts given in the RC3 catalog (z = 0. 014300, 0.024307, and 0.018 126, respectively: Ide Vaucouleurs et alJl99lTl : IRAS 23128- 
5919: luminosity distance computed bv lCharmandaris et a,lJ (I2002T) . All distances assume Ho = 75 kms -1 Mpc -1 and qo = 0.5 
where relevant. 

6 Morphological types are from the RC3 catalogue (Ide Vaucouleurs et al.lfl99 lT) . with additional indications for interacting 
("int.") or merging system. 

c Nuclear types g ive an indication of the nuclear activity inferred from optical spectra. The types for the normal spiral galaxies 
are as listed by I -loussel et alJ i2001bl : see their Table 1 for references). For a subset of the starburst s ystem s, references 
arelKewlev et all ll200lt NGC 253. N1808, IRAS 23128-59191. IStanfordT Il99ll : NGC 520), j^ejUeux et alJ Jfl995t NGC6240), 
llGrri^t^dTi^9^ : Arp 220). The types for the remaining starbursts are those listed i n the NED database . 

d Global infrared (8 — 1000 /im) luminosity computed from the IRAS fluxes following ISanders fc Mirabell <ll996f) . 

e Data available from ISOCAM observations and used in our analysis: "LW" for data obtained in the broad band filters LW2 
and LW3 (5.0 — 8.5 /im and 12 — 18 /^m), "CVF" for 5 — 17 /im spectrophotometric imaging at R ~ 40. 

' Observing program. SLglx: P.I. G. Helou; CAMbarre: P.I. C. Bonoli; CAMspir: P.I. L. Vigroux; Virgo: P.I. J. Lequeux; 
Zzcvfcam: P.I. D. Cesarsky; CAMACTIV: P.I. Mirabel; IM SPJ5BG: P.I. R. M aiolino, 

9 Refere nces for initial publicatio n of the ISOCAM data: (1 ) iDale et alJ (l2000l): f2llRoussel et"ai1 (l2001al): OIlRoussel et alJ 
fcOOlbl): r4)lBoselli et alJ (ll998Tl: (5) ISauvaee et all (Il996r) : (6 ) iMalhotra et alJ Jl996l ): (7) IForster Schreiber et alJ (l2003li T 
f8l lLaurent et alJ i200oK 7 (^ ICharmandaris et alJ J2002ft : (10) iTran et all fcOOll) 

t This galaxy also is a LIRG, but about half its total mid-infrared emission arises outside the central regions selected here. 



been published anywhere else; they are briefly presented 
in Appendix iBl 

The sample can be divided in four parts, in order of 
increasing star formation activity: 

— disks of spiral galaxies, which form stars in a quies- 
cent fashion and for which the relationship between 
the MIR emission and the SFR, as der ived from Hcv 
line m easurements, was discussed by iRoussel et all 
feOOld) : 



— the more active circumnuclear regions of a subsamplc 
of the same galaxies, for which it was possible to es- 
timate an extinction-free SFR from a combination of 
near- infrared H recombination lines and Ha; 

— nearby starburs t galaxies, three of which were stud- 
ied in detail by IForster Schreiber et all l|200ll |2003|) 
in their dust and fine-structure line emission; 

— luminous and ultraluminous infrared galaxies (LIRGs 
and ULIRGs, with 10 11 < L m < 10 12 L and L m > 
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10 12 L qj respectively 1 ) take n from the samples stud- 
" tOOd) . iTran et alJ tOO^i and 



ied bv ^Laurent et al 
ICharmandaris et al.l 



Most objects are purely star-forming systems; a few 
harbour an active galactic nucleus (AGN) revealed by op- 
tical spectroscopy, which is however known or suspected 
to contribute but negligibly to the MIR emission, as well 
as to the hydrogen recombination lines within our pho- 
tometric apertures, which are very large with respect to 
the angular size of the nucleus. Detailed notes and dis- 
cussion of individual cases are given in Appendix The 
LIRGs and ULIRGs will also be referred to as starbursts 
throughout the paper. Tabled lists the non-disk objects 
(circumnuclear regions of spirals and starburst systems). 
The table gives some general properties, the observation 
mode, and the original observing program to which they 
be long. Similar details for the spiral disk sample are given 
bv lRoussel et alJ <|2001albl) . 

In spiral galaxies, the low brightness disks and central 
regions are distinguished by different ratios of flux den- 
sity in the LW3 and LW2 filters, or /12-I8 »m/ f5-8.5»m 
color: the disks typically have ratios of ss 1 while the 
circumnuclear regions usually exhibit a color excess si; 
naling more active s tar formation 2 (e.g. iDale et al1l200 
iRoussel et al.ll2001bh. We adopte d the m easurements for 
Li: 



disks reported bv IRoussel et all l)2001ch . Briefly, these 
were obtained from the integrated MIR and Ha fluxes 
by subtracting the contribution from a core region and 
accounting for flux dilutio n effects of the ISOC AM point 
spread function (PSF; see [Roussel et al.ll20f)la|) . The size 
of the excluded area was dictated by the Ha data exist- 
ing in the literature (fluxes in given apertures, or maps). 
In the few cases where the Ha aperture is smaller than 
-Ccnr, the size of the circumnuclear regions fitted on MIR 
brightness profiles, it was ensured that the resulting disk 
/i2-i8A 1 m//5-8.5^m color was close to unity. 

Altogether, the sample covers more than five orders 
of magnitude in Lyman continuum photon flux density 
^Lyc (Sect. 0J. The latter implies five orders of magni- 
tude in quasi-instantaneous SFR surface density, assum- 
ing that the same stellar initia l mass function applies to all 
objects (e.g. lKennicuttlll998|) . Our sample is admittedly 
not complete in any sense and is restricted to near-solar 
metallicities. For our purposes, it should however provide 
a sufficiently representative ensemble since the primary 
samples were constructed with different criteria and aims. 
Since star-forming systems, both Galactic and extragalac- 
tic, have remarkably similar MIR spectral energy distribu- 
tions (SED) in terms of broad features and continua, we 
are confident that we are not introducing any bias by se- 
lecting particular galaxies. The sample was only shaped by 



1 Lir, = I/8-1000 nm is the total infrared luminosity 
computed from IRAS f luxes following the prescription of 
ISanders fc Mirabell <1996ft . 

2 More generally, such a color excess can also be due 
to an AGN heatin g a surrounding nuclear dust torus (e.g. 
lLaurent et aIl l2uOCn . 
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Fig. 1. Bandpasses of the MIR indicators of star forma- 
tion activity. The transmission profiles of the ISOCAM 
LW2 and LW3 filters and the narrow band used for the 
monochromatic 15/im,ct continuum measurements are 
shown by the grey lines. To illustrate the typical spec- 
tral features covered by these bandpasses, the spectra of 
the starburst galaxy M 82 obtained at R ~ 500 - 1000 
with the ISO SWS and at R ~ 40 with ISOCAM within 
the SWS field of view are plotted in the top and bottom 
pane ls, respectively (from iForster Schreiber et al.l l200ll 
120031) . Identifications of the emission features are labeled 
on the SWS spectrum. 

the availability of adequate data. Further details including 
notes on individual sources are given in Appendix lAl 

3. Star formation diagnostics 

Tables|21and|21report the data for circumnuclear regions of 
spirals and starburst systems that we used i n our analysis. 
The da ta for spiral disks are described bv iRoussel et alJ 
fcOOld) . We reduced and analysed the MIR maps and 
spectra of the whole sample in a homogeneous way. 

3.1. Mid-infrared 5 — 8.5 fim and 15 [xm emission 

The shape of the 5 — 11 /im SED is observed to be 
nearly invariant in star-forming galaxies and in a vari- 
ety of Galactic sources while at A ^ 11 fim, the sub- 
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Table 2. Dust emission of circumnuclear regions of spiral galaxies and starbursts 



5 



Source 


Region a 


Area" 


f b 

75—8.5 (im 


fl2 — 18 


b 


fl5 /j. 


C 








T W9 


V v r 


LW3 


CVF 


LW 


' CVF 




(arcsec) 


(PC 2 ) 


( mJ y) 


frr, l^A 

(mjyj 


(mjv) 


(mjv) 


(mjv) 


(mjy) 


Circumnuclear regions of spiral galaxies 














IN vi v yoo 


9 1 } £ 18 9 


^ 97 v 1 n" 

O.Z I A 1U 


^9/1 

OZ^i 




673 




628 




at r^r 1 1 no7 
IN Illy f 


4-J 


7 qg w 1 n^ 


1 98n 


(A /I QQ^ 
^14oy ) 


1692 


(2200) 




1641 


IN loDO 


/in n Q7 i 
4U.U, ol.L 


7 ofi v i n^ 

/ .ZD X 1U 


1 Q79 

iy / z 


(A 5Z$i7\ 


3103 


(3451) 




2940 


Mrr^ /i i no 

IN ^jtU 41UZ 


on /i oq 7 


4.oy x iu 


> 571 




> 1628 




1673 




\rpp /j ono 
IN <±Zt)o 


qo /l fin 
OZ.4, O.U 


/i *3Q v 1 n^ 
4.oy x iu 


7fi 




147 




134 




IN 40y 1 


cfi c on /i 

00. 0, ZU.4 


77 v 

Z. it X IU 


^49 




759 




582 




wnp ^1 Qd 


Qfl 


G87 y If) 6 
O.O I A IU 


1 88^ 


^zuyu ) 


2021 


(2292) 




1382 


i,l VI V. / UZiOU 


41.2 38.4 


01 x 1 5 


z ( oo ) 


^1 OQ 


(> 3589) 


5468 




4369 


IN oy4o 


qc q qo n 
oO.O, OZ.U 


D.l x X IU 




1 ^0^ 


(> 1885) 


2174 




1551 


ivrr 7^9 


Ol C 1 Q Q 


o 7Q Y in 6 

z . ( y a _lu 


1248 




> 2316 




2067 




ivrr 7771 


99 A 90 ft 

ZZ.4:, ZU.U 


O Cjfl y 111 7 
Z.OU A _LU 


314 




387 




270 




Sffirhnrtif nnl (TV) p i 


















NGC 253 


15.9, 14.0 


2.27 x 10 4 




yJOOO 




22 755 




23 913 


NGC520 


17.7, 15.3 


3.36 x 10 6 


^78 




737 


(670) 




560 


NGC 1808 


25.9, 24.4 


1.28 x 10 6 


2645 


(2450) 


4077 


(4448) 




3426 


NGC 3034 


30.0, 28.6 


1.65 x 10 5 




25 812 




62 366 




60 000 


IC342 


17 x 17 


7.40 x 10 4 




1404 




3527 




3035 


LIRGs/ULIRGs 


















NGC 3256 


20.0 


1.03 x 10 7 




1543 




2960 




2634 


NGC 6240 


total (3.) 


1.57 x 10 6 


190 


(235) 


767 


(784) 




801 


IRAS 23128-5919 


total (4.2) 


1.08 x 10 7 


120 


(119) 


331 


(468) 


337 


(443) 


Arp 220 


total (2.) 


2.11 x 10 6 


191 


(196) 


765 


(869) 




1170 



a "Region" refers to the rectangular dimensions or diameter of the photometric aperture; the first value is the aperture, and 
the second value is the deconvolved size, which was used for the area normalization as given in "Area". 

b Broad-band flux densities from ISOCAM observations measured either through the LW2 and LW3 filters or computed from 
CVF spectra accounting for the filter transmission profiles. Lower limits for LW measurements indicate that the nucleus is 
saturated in the maps; in which case measurements from the CVF spectrum were preferred. Values unused in the present 
analysis are enclosed in parentheses. From unpublished previous analysis, photometric errors are always dominated by 
incompletely corrected memory effects in the LW2 filter, and in the LW3 filter for sources brighter than w 200 — 500 mjy. 
Taking these memory effects into account, it is estimated that average errors are 10% and 2 0% in LW2 and LW3 , respectively; 
individual errors for relatively bright galaxies may be as high as 20% and 30%, respectively jRoussel et al.l2001al) . In addition, 
flux calibration uncertainties are of the order of 5% (ISOCAM handbook). 

c Monochromatic 15 /im continuum flux density derived from LW2 and LW3 data bas ed on Fig. E3 as explain ed in Sect. 13.11 
or measured directly from CVF spectra. Calibration errors for CVFs are discussed bv lBiviano et all (ll998aT) . 



stantial drop in the most quiescent galaxies, with con- 
tributions by minor aromatic features, contrasts with 
the increasingly strong and steep continuum of VSGs 
in more active sources (see Tielens 1999 for a review; 
see also e.g. Bou l anger et al. 1998: iHelou et al.l J2 000: 
Uchida et al.l2000HLaurent et al.l2000HSturm et al|2 000: 



3.6 /im PAH complex where the underlying continuum 



Roussel et all l200ld iForster Schreiber et all [2003). In 
ULIRGS, extinction effects can be large enough to distort 
the shape of the PAH complexes especially by the suppres - 
sion of the 8.6 /xm feature (e.g. iRigopoulou et alj|l9 99). 
PAH bands are generally very weak or absent in spectra 
of pure Hll regions and AGNs, a fact usually attributed 
to the destruction of band carriers in hard and intense ra- 
diatio n fields (e.g. IRigopoulou et al1ll999t lLaurent et al.1 
2000;, and references therein). 

We focussed on two bandpasses sampling as indepen- 
dently as possible the PAH and VSG emission. The LW2 
band (5 — 8.5 ^m) encompasses the prominent 6.2, 7.7, 



stellar activity (e.g. 


Rigopoulou et alJll99Qt lLaurent et alJ 


2000: iLu et al.ll2003 


). To probe the VSG emission, we pre- 



ferred to define a narrow interval measuring the monochro- 
matic flux at 15 /Ltm rather than use broadband measure- 
ments through the LW3 filter. The LW3 bandpass in- 
cludes the strong 12.7 /im PAH as well as minor features 
at 13.55, 14.25, and 15.7 //m that pro bably dominate the 
SEP at low star f ormation levels (e.g. [Sturm et aliboOOt 
iHonv et all l200ll iRoussel et all 12001^ . We defined the 
15/Ltm narrow band as a top- hat profile filter with unit 
transmission between 14.8 and 15.2 /Ltm, maximizing the 
VSG contribution by avoiding known PAH features and 
other possible emission lines. We note that line emission 
is not expected to contribute significantly to LW2 and 
LW3 measurements. For instance, spectra of the ISO Short 
Wavelength Spectrometer (SWS; Ide Graauw et al.l ll996) 
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Table 3. Line fluxes and extinction for circumnuclear regions of spiral galaxies and starbursts 



Source 


Aperture a 


Fx 


b 






A v c 


log[QLyc] d 


References e 




(arcsec) 


(HT 17 Wm" 


2 ) 


(10- 16 Wm" 2 ) 


(mag) 


(logls- 1 ]) 




Circumnuclear regions of spiral galaxies 
















IN yoo 


23.0, 4X4 








10.2 


2.5 


53.57 


1, 2 


in v < v ±uy ( 


1 r » 








39.2 


1.9 


53.55 


3, 4 


X\ vi v 1O00 


/inn 9^ ^ 


Rr'V • 

di y . 


19.6 




35.8, 28.4 


2.9 


53.95 


5, 6 


Kir^r 1 a 1 no 


starburst 


r Ap J r>r 7 


12.1/6.0 




7.2 


53.55 


7 


IN vjL> QZito 


Q9 /I 




10.7 




0.76 


3.8 


52.59 


8, 9 


. n v i \ ^toy i 


^fi ^ 9*3 ^ 


or 7 . 


3.76 




8.59, 7.18 


2.4 


53.45 


10, 11 


IN V < v iJ-Ly'i 


on 








38.6 


2.8 


53.29 


12, 13 


NGC 5236 


41.2, 23.5 


Br 7 : 


19.6 




75.7, 60.6 


1.7 


52.79 


5, 14 


NGC 6946 


35.3, 23.5 


Br 7 


6.1 




8.64, 5.80 


3.5 


52.58 


5, 15 


NGC 7552 


91 8 14 x 9(1 


Bra : 


44.0 




30.3, 29.4 


2.1 


53.75 


1, 16 


NHP 7771 

IN V I V., / / / J. 


o Ldi yj u.i bl i nit; 


Pa/3/Br 7 


: 4.0/1.5 




5.2 


53.91 


17 


S //i T'hii t"c/ rtnl t~t htq & q 


















NGC 253 


15, 2.4 x 12 


Br 7 : 


91.6 






8.5 


53.13 


18 


NGC 520 


17.7, 6x8 


Br 7 


1.9 




0.52, 0.15 


7.4 


54.03 


19, 20 


NGC 1808 


20 (total) 


Br 7 : 


31.0 






4 


53.72 


21 


NGC 3034 


30.0 










52 (MIX) f 


54.09 


22 


IC342 


17 x 17, 14 x 20 


Br 7 : 17.0; Bra 


70.0 




5.2+ 


52.48 


16, 23 


LIRGs/ULIRGs 
NGC 3256 ? 


a 20 (total), 3.5" x 3.5" 


Br 7 : 


15.0 






5.3 


54.54 


24, 25 


NGC 6240 


total 


Br 7 : 3.1; Pa/3 


4.1 




10.1 


54.91 


26, 27 


IRAS 23128-5919 


total, S. nucleus 


Paa 


8.0 




1.76, 0.93 


3.0 


54.74 


28, 29 


Arp 220 


total 


Bra : 21.0; Br 7 


0.59 




40. L 


55.33 


30, 31 



a Rectangular dimensions or diameter of the photometric aperture. When two values are given, the first one refers to the 
primary flux (Ha when present) and the second one to the hydrogen line decrement used to estimate the extinction. The 
aperture adopted for the size normalization is given in Table [5] 

b Observed line fluxes used to derive the extinction and/or the ionizing photon flux. When two Ha fluxes are given, they 
correspond to the two apertures used for photometry and extinction estimation, respectively. 

c Derived or adopted extinction for a uniform foreground screen (UFS) model except when "MIX" indicates that a homoge- 
neous mixture of dust and sources is assumed. 
Derived or adopted intrinsic Lyman continuum photon rate. 

e References for line flu xes and extinction values (f urther details are given on the derivation of the data in append ix lAl : 
(1) Ha+[N II] map from lHameed fc Devereip] ljl99flfr (in NED); (2) Extinction f rom H/3/Ha oflVeron-Cettv fc Veronl l|l98fih 
in 4" x 4"; (3) Ha+[Nll] map graciously provided by T. Storchi-Bergmann dStorchi-Bergman et al.l|l996| ); (4) Average 
extinction along the starburst ring derived from the data of iKotilainen et all (l200Cl) ; (5 ) Br 7 flux from|Pjr xley_et_alJ il988l) 
in an effective aperture of 23. 5"; (6) Ha+[Nlll ma p graciously provided by M. Naslund (|Kristen et alJ ll997l : lLindbradll 999): 

i 7) Br 7 and Pa/3 fluxes from lRoussel et all l|200 3|) in ~ 5" x 5" (m ost of starburst enclosed); (8) Paa map from l Boker et all 
19991) (in NED); (9) Hq+[Nul map from lKoopmann et"aH fcOOll) (in NED); (10) Br 7 flu x fromlPuxlev et alJ Jl990ll l in an 
effective aperture of 23.5"; (11) Ha+[Nll1 ma p graciously provided by A. Garcfa-Barreto §Garcfa-Barreto^^ilJ^)99) ; (12) 
Ha+[Nll] map from iGreenawalt et al] {l998); (13) Ave rage extinction tow ards Hn regions from pcoville et alj (|200l| ); (14) 
Ha map graciously provid ed by S. Ryde r via A. Vogler jRvder et aljll995l) : (15) Ha map fr omlLarsen fc Richtlerl dl999T) (in 
NED); (16) Bra flux from lVerma et al] J2003D in 14" x 20"; (17) B r 7 and Pa/3 fluxes fromlDale et al] (120031) in s a 5" x 5" 
(most of starburst ring enclosed); (18) Br 7 fl ux and extinction from [Engehb racht et all dl99a) : (19) Bro- flux from IStanfordl 
19911) in 6" x 8"; ( 20) Ha+[Nll] map from lHibbard fc van Gorkoml (I1996T) ( in NE D); (21) Br 7 flux and extinction from 



|Krabbee^d 



,1994); (22) Extinction and QLyc from Forster Schreiber et al. (2001); (23) Br 7 map provided by T. Boker 
19971) : (24) Total Br 7 flux from iMoorwood fc Olival lll994l) : (25) Extincti on from Br 7 /Pa/3 of iDovon et al] 



Boker et al] | ... _ . _ 

19941) in 3.5" x 3.5"; (26) Br 7 flux from lRieke et al] lll985l) in 8.7"; (27) P a/3 flux from lSimpson et al] ]l996l) in a slit of 



width 1.5"; (28) Paa flux of the southern nucleus from [ Kawar a et ap 1119871) : (29) Ha fluxes of both nuclei form |Duc_et_al] 
Jl997f) in a slit of width 1.3"; (30) Bra flux from lSturm et al] <ll996j) in 14" x 20"; (31) Br 7 flux of lGoldader et al] (11995 
in a slit of width 0.75". 

t Values are for the GC extinction law of iLutzl (ll999aT) when we made use of the Bra line at 4.05 /jm (see Sect. 13.21 . 



at R ~ 1000 show that the strongest lines falling within 
the LW2 and LW3 bandpasses for M82 and NGC 253 
are [Aril] 6.99 /xm [Nell] 12.81 /xm, and [Neill] 1 5.56 
l|Sturm et alJl200a lF5rster Schreiber et al]l200l|) ; we de- 



termined that they account for only « 1% and 3% of the 
LW2 and LW3 flux densities, respectively, in both galax- 
ies. 
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•/ N6946 


0.5 






0.0 







12-lB/im 



5-8.5/^m 



Fig. 2. Empirical relationship used to estimate the 
monochromatic 15 (im continuum flux density from broad- 
band LW2 and LW3 measurements. The black circles rep- 
resent circumnuclear regions of spiral galaxies and star- 
burst systems for which both broadband and CVF data 
are available. Diamonds represent galaxies with only CVF 
data, for which broadband fluxes were synthesized from 
the spectrum. The line shows the least-squares fit obtained 
as explained in the text, y = a — b/x, where x and y are 
the abscissa and the ordinate. The error bar indicates the 
mean and ler dispersion in /i2-i8 ( um//5-8.5^m color of the 
sample of spiral disks, and the /i5 Mm ,ct//i2-i8/*m ratio 
and measurement uncertainty for the disk of NGC 5236. 
The disk of this galaxy was chosen because its spectrum 
has the highest signal-to-noise ratio of the sample. 



For NGC 1097, this is probably because our correction 
of the instrument's memory effects is most uncertain at 
A > 1 6 ^m. The correction algorithm I Coulais fc Abergel 
200fl| requires knowledge of the detector's response dur- 
ing previous exposures. Since CVF spectral scans were 
performed in order of decreasing wavelength, and the illu- 
mination history of the detector prior to our observation 
is unknown, the spectrum is most affected at long wave- 
lengths. The CVF spectrum of IRAS23128-5919 has the 
lowest signal to noise ratio in the sample and is therefore 
less reliable; we do not use it in what follows. 

Six galaxies of our sample were observed only in broad- 
band photometric mode, so we derived an empirical con- 
version between /12-18/mi an d /i5^m, ct as follows. We re- 
lied on the interpretation that the emission in the LW2 
band is dominated by PAHs, and that the emission in the 
LW3 band is produced mainly both by PAHs (or akin par- 
ticles) and by VSGs whose flux variations behave differ- 
ently (see Sect.GJ. Under this assumption, /i5^ m)C tj which 
is covered by the LW3 filter, should also contain contri- 
butions from these two species, albeit in different propor- 
tions, and we should expect it to be related to the broad- 
band measurements by a simple function. Figure [5] shows 
the exact values of this function taken by the objects for 
which /5-8.5/jm and fi2-isnm were measured from LW2 
and LW3 maps, and /ls^m.ct from CVF data. Assuming 
that the emission covered by the LW3 filter is the sum 
of a component scaling linearly with the emission bands 
covered by LW2 and of a second component scaling lin- 
early with the emission measured by /15 Mmi c t , it is easy 
to show that /ls^m, ct//i2-i8^m can be represented by an 
affinc function of l/(/i2-i8^m//5 — 8.5 /im ) • 



Throughout this paper, we refer to the narrow 15 /im 
bandpass as "15 /im, ct" and adopt the notations "5 — 
8.5 /im" and "12-18 /im" for LW2 and LW3. The effective 
bandwidths are 0.53, 16.18, and 6.75 THz, respectively. 
Figure ^ shows the corresponding wavelength ranges and 
transmission profiles on the SWS spectrum of M 82 and 
on the lower resolution ISOCA M spectrum extracted 
within the SWS field of view l|F6rster Schreiber et alJ 
l200lM2003h . 

We obtained the /5-8.5/im flux densities directly from 
LW2 observations when available, or computed them from 
CVF spectra accounting for the LW2 filter transmission 
profile. Values of /5-8.5/im derived from LW2 and CVF 
data for all the sources observed in both modes agree 
within 20%; differences may be attributed in part to pos- 
sible residuals f rom ghosts, flat field, and straylight i n 
the CVF data l|Biviano et al.l Il998all4 lOkumural^OOol) . 
Extrapolation of the spectra between w 16 — 17 /im and 
18 /im is necessary to compare the /i2-i8pm fluxes derived 
from LW3 and CVF data; however, for all the galaxies ob- 
served in both modes, we find differences of less than 20%, 
except for NGC 1097 and IRAS23128-5919, whose CVF 
data overestimate /i2-i8£im by 30% and 41%, respectively. 



Excluding Arp220, whose CVF spectrum suggests se- 
vere extinction effects, with high opacity from amorphou s 
silicates at 9.7/mi and 18 /im l|RigODOulou et al] |l999). 
and NGC 1097 whose CVF spectrum is affected by resid- 
ual memory effects, we used the least-squares fit to these 
data to assign a /i5 Mm . ct//i2-i8/mi ratio to the circumnu- 
clear regions of spirals and starbursts without CVF data. 
Figure shows also data for galaxies observed only in 
CVF mode, whose broadband fluxes were simulated from 
the spectrum (diamond symbols). Except for Arp220 and 
NGC 1097, all data points are within 15% of the fitted re- 
lation. For spiral disks, we applied a uniform conversion 
j ustified by their sm all dispe rsion in fi?-is„rn/fR-s s».m 
I Roussel et al.lho01(4 see alsolDale et al Enrol l200ll We 
used the ratio fi5^ m , ct//i2-i8/im = 0.38 measured on the 
disk of NGC 5236 which has the best quality CVF spec- 
trum among the disk sample. 

We did not correct MIR flux densities for extinc- 
tion. Relative to the optical V band (5500 A), the ex- 
tinction in magnitudes is ver y small, with A^-h 5 , UTTI I Ay 
and A 15fim /A v < 0.06 (e.g. lDrainelll989t lLutz|]l999a|) . 
Extinction effects on the relationships studied in this work 
will be discussed in Sect. 0] 
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3.2. Star formation rate indicators: H recombination 
lines 

To estimate SFRs, we used H recombination lines collected 
from the literature, which provide primary diagnostics and 
allowed us to derive the nebular extinction and correct for 
it. We converted all fluxes to a common reference quan- 
tity, the production rate of Lyman continuum photons 
QLyc- We took care that consistent apertures were used to 
measure the dust and hydrogen line fluxes. Limitations of 
available data and the assumptions on physical conditions 
made in deriving QLyc inevitably lead to appreciable un- 
certainties. We emphasize however that uncertainties for 
individual sources of even a factor two affect but little our 
conclusions, as will be discussed in Sect. 0] 

For spiral disks, we used the total and circumnuclcar 
Ha + [Nil] AA6548 ,6583 A fluxes correc ted for Galactic 
extinction listed bv lRoussel et al . (2001c) to obtain disk- 
only fluxes. W e then derived in trinsic Ha fluxes following 
the precepts of lKennicuttl l|l983f) applicable to H II regions 
in spiral disks, correcting for an average 25% contribu- 
tion by the [N n] lines and an average inter nal extinction 
Aun = 1-1 mag (see also the discussion bv lRoussel et all 
200l]J). 

For the circumnuclcar regions of spirals with two or 
more H line measurements, we derived the extinction by 
least-squares fits to the ra tios of observed fluxes to in trin- 
sic line emissivities from IStorev fc HurrJ ill 995ft . We 
assigned equal weight to the ratios given the difficulty of 
determining the uncertainties for the inhomogeneous col- 
lection of H line data. We obtained the resulting Ql Y c by 
averaging the individual Qhyc values derived from each 
dereddened lin e flux, taking the total H recombination 
coefficient from lStorev fc Hummerl l|l995lh In some cases, 
relevant line data were not available or too uncertain, so 
we adopted published values of extinction insofar as deter- 
mined from H lines. Since the extinction was sometimes 
derived in a region smaller than our photometric aper- 
ture, and the assumption of uniform extinction through- 
out kiloparsec scales is probably wrong, extinction cor- 
rections may introduce a non-negligible dispersion in the 
relations shown in Sect. 21 

Whenever circumnuclcar Ha fluxes included the satel- 
lite [Nil] lines, we applied the same correction factor 
of 0.75 as for the disks, the validity of which we veri- 
fied as much as possible based on publis hed spectroscopy 
from v a rious sources. The com pilations of | K ennicut t et alJ 
l|l989j) . iKennicuttl <|l992j) and Ijansen et all l|200f| shew 
that while there is significant overlap in [N n] A6583 A/Ha 
ratios between disk H n regions and nuclei of spiral galax- 
ies, many (~ 50%) can exhibit much higher ratios, which 
can be explained by shock heating (e.g. by supernova rem- 
nants) or by a non-thermal Liner/Seyfert contribution. 
This line ratio increase is however observed in the im- 
mediate vicinity of nuclei; our circumnuclear regions are 
generally much larger so that this effect is not expected 
to be important. 



W e employed the ex tinction law of ICardelli et alJ 
l)l989h at A < 3 ^m and of lLutd l)l999ah at A > 3 /^m. We 
adopted a uniform foreground screen model (UFS) for the 
geometry of the sources and obscuring dust. The limited 
number of H lines considered for each galaxy prevented us 
from constraining the extinction model, but computations 
for a homogeneous mixture of dust and sources ("MIX" 
model) imply QLyc values differing by at most 56% (on 
average 18%) from those of the UFS model 3 . We assumed 
that the Hn regions in all sample galaxies are ionization 
bounded, optically thick in the Lyman lines and optically 
thin in all others (case B recombination), and adopted 
electron density and temperature of n c — 100 cm~ 3 and 
T c = 5000 K. These n c and T c were found representa- 
tive for a sample of starbursts observed with SWS by 
iThornlev et al.1 ifeoOO'l. including most starbursts in our 
own ISOCAM sample. Higher values up to n c = 10 4 cm~ 3 
and T c = 10 4 K may be more appropriate for H n re- 
gions in_jwnnal spiral galaxies of near-solar metallicit y 
(e.g. ISmithllT975t IShaver et aljFflsi ICiveon et ahlEflol. 
However, the computations of IStorev fc Hummerl l|l995h 
imply variations of the relative emissivities of 13% on 
average (25% at most) for the lines considered here, lit- 
tle affecting the extinction estimates (< 0.5 mag for the 
UFS model). The mean increase in the derived Qi,y C be- 
tween n c = 100 cm~ 3 , T c = 5000 K and n c = 10 4 cm" 3 , 
T c = 10 K is 36% (maximum 55%), mainly driven by the 
variations of the total H recombination coeffi cient qr with 
Tf. (a R depends only weakly on ro e ; IStorev fc Hummerl 
1995). In addition, we do not consider individual bright 
Hn regions, but the total emission from large areas en- 
compassing many star formation complexes. The average 
values of n c and T e are thus expected to be much lower 
than in resolved Hn regions. 

3.3. Size normalization 

To obtain scale- and distance-independent quantities, we 
normalized each measurement by the projected surface 
area and expressed the results in L Q pc -2 (denoted here- 
after £5-8.5 Aim) Si5/im,ct, ^Lyc)- We chose these units 
specifically to avoid artificial correlation due to scale ef- 
fects whereby the brighter (larger) galaxies tend to be 
brighter at all wavelengths. Normalizing all three quanti- 
ties by the surface area eliminates dispersion from uncer- 
tainties in distance estimates. The QLyc values were trans- 
formed into Lyman continuum luminosities ^Lyc assuming 
an average ionizing photon energy of 16 eV. We empha- 
size that although these quantities are formally equivalent 
to surface brightnesses, they are not intended as such and 
will be referred to as "size-normalized luminosities." 



3 The observed and intrinsic line fluxes Fx and F® are related 
through Fx = F° e" r * for the UFS model and Fx = Ff [(1 - 
e~ Tx )/rx\ for the MIX model, where tx is the optical depth 
of the obscuring material and the corresponding extinction in 
magnitudes is Ax = 1.086 Tx- 
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We normalized fluxes of spiral disks by the circu- 
lar area of diameter -D^, the major axis length of 
the B-band i sophote /xr = 25 magarcs ec -2 (from the 
RC3 catalog; Ide Vaucouleurs et al.lll99lj) . This area en- 
closes all detected MIR emission as defined by the LW2 
isophote at 5^iJyarcsec -2 , the typical depth reached 
in these ISOCAM data, with isophote diameter ratios 
£>5~j y 5Mm /D^ in the range 0.35 — 1 , depending on the ga s 
richness and inclination of galaxies l|R.oussel et alJEoOla^ . 
For galaxies with no available Ha map, it was verified 
that the aperture of the integrated Ha flux is larger than 
or comparable to the spatial extent of the MIR emission. 
The mismatch between the optical diameter and the ac- 
tual sizes of the MIR and Ha emitting regions will in- 
troduce either a "correlation bias," whereby points move 
along a line of slope 1, or scatter in our relationships, de- 
pending on how well the MIR and Ha emission trace each 
other and are covered by the photometric apertures. 

For the circumnuclear regions and starburst cores, we 
have striven to adopt as well a uniform and well-defined 
surface quantity. We extracted the size of the MIR source 
from azimuthally-averaged surface brightness profiles fit- 
ted with a gaussian, after bringing all images to a common 
gaussian PSF (see Appendix lA"|l . We chose a photometric 
aperture of 2.5 times the fitted half-power beam width 
(HPBW) in order to measure the total flux to a good 
aproximation; after deconvolution by the size of the PSF, 
the aperture is also used for area normalization. For the 
sources whose geometric structure is ill-suited to such a 
definition, we made exceptions explained in Appendix lAl 
In particular, for the most compact sources (with re- 
spect to the ISOCAM angular resolution), NGC4293, 
NGC6240, IRAS23128-5919 and Arp220, we obtain in 
this way intrinsic sizes that are comparable to or smaller 
than the angular resolution, and therefore highly uncer- 
tain. Hence, for these galaxies, we did not use the fitted 
MIR sizes for area normalization, but starburst sizes de- 
rived from high-resolution observations in the literature. 



4. Results 

4.1. Calibration of MIR dust emission as star 
formation diagnostic 

Figure 01 shows the relationships between the observed 
^i5^m,ct and £7 ^ and the derived SL yc for our sam- 
ple galaxies. The immediate result is that both MIR 
quantities constitute reliable SFR tracers over many or- 
ders of magnitude in SLyc- Thus, in a general sense, 
our results extend those previo usly found for spi- 
ral disks by iR.oussel et all l|2001ch and for NGC 253, 
NGC1808, and M82, i nclud ing individual regions, by 
iForster Schreiber et al Our enlarged sample re- 

veals however interesting differences in the behaviour of 
the MIR tracers and between various source classes. 

In the Si5 Mm ct — J^Lyc diagram, the monochromatic 
15 /im continuum emission is directly proportional to the 
ionizing photon luminosity within the error bar on the 
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Fig. 3. Empirical relationships between the MIR dust 
emission and the Lyman continuum luminosity as a mea- 
sure of the star formation rate. The quantities plotted are 
size-normalized luminosities, with different symbols used 
for the spiral disks (crosses), and the circumnuclear re- 
gions of spirals and starburst systems (diamonds), (a) 
In the log(£i5 Mm , ct ) — log(SLyc) diagram, the solid line 
indicates the power-law least-squares fit performed on 
the circumnuclear regions and starbursts only, excluding 
the disks, NGC 5194 and Arp220. The short-dashed line 
shows the fit result where the index was fixed to unity, and 
the long-dashed line shows a similar fit for disks alone, (b) 
In the log(S7 Aim )— log(SLyc) diagram, the conventions are 
the same and fits were performed including the disks, but 
excluding NGC 253, NGC 3034 (= M82), NGC 6240 and 
Arp 220. 



power-law index, but with two different normalizations 
corresponding to two distinct regimes, each spanning 2 — 3 
orders of magnitude in SLyc^ 1) quiescent spiral disks, 
with low Lyman continuum luminosities per unit pro- 
jected area El yc 10 2 L Q pc~ 2 ; and 2) moderately to ac- 
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tivcly star-forming regions in the central ^ 1 kpc of spiral 
and starburst galaxies, with an activity level characterized 
by ^Ly C iS 10 Lq pc , but perhaps excluding extreme 
environments with ^Lyc ~ 10 5 — 10 6 L Q pc~ 2 . Although 
the transition between the two regimes is certainly grad- 
ual, our data set includes but one object sampling it. It 
occurs approximately at the level of star formation activ- 
ity seen in the inner « 90" plateau of M51, which has 
a color /i2-i8/im//5-8.5/jm ~ 1-1- The offset between the 
normalization of disks and that of starbursts is nearly a 
factor 5, and cannot be caused by the different methods 
applied to estimate ionizing photon fluxes. Assuming di- 
rect proportionality in each separate regime, we obtain: 

log(Si 5 „ m , ct ) = log(S Lyc ) - 2.28, log(S Lyc ) < 2, (1) 
) = log(E Lyc ) - 1.60, 2 < log(E Lyc ) < 5, (2) 

where all size- normalized luminosities are in Lq pc -2 . The 
fits are shown as dashed lines in Fig. [3] The dispersions 
are, respectively, 0.22 dex for the disks (factor w 1.6), and 
0.18 dex for the galactic centers and starburst systems 
(factor 1.5), from which we have excluded NGC5194 and 
Arp220 at the extremes of the I^Lyc range. Linear least- 
squares fits, where the power-law index is let as a free 
parameter, yield an exponent of 1.01 ±0.07 for disks and 
0.97 ± 0.06 for starbursts. 

The break between spiral disks and more active regions 
can be interpreted easily. In disks, the density of the radi- 
ation heating the dust is too low for the VSG continuum 
to be significant at 15 /xm. In such conditions, the contin- 
uum starts rising at longer wavelengths where larger dust 
grains at lower temperatures re-emit the energy absorbed 
from the relatively diffuse radiation field. The emission 
detected at 15 /im is then dominated by PAHs or a re- 
lated family of particles. Above a certain threshold in ion- 
izing radiation density, VSG heating becomes more effi- 
cient such that the continuum starts to make a significant 
contribution at 15 /im, and another regime prevails. The 
break thus signals the onset of VSG emission at sufficient 
star formation densities. Since linear correlations repre- 
sent adequately the data in both regimes, we infer that 
the respective contributions from each dust species vary 
almost linearly with Lyman continuum luminosity. The 
data of Arp 220 suggest a flattening of the relationship at 
the most extreme densities; its Si5 Mm!Ct lies about a factor 
of 4 below a simple extrapolation of the linear correlation 
seen in the other starbursts and the circumnuclear regions. 
Since the MIR spe ctrum of Arp 220 shows some signs of 
high optical depth ijRigopoulou et al]ll999|) . this damping 
of the dust emission can easily be caused by extinction 
effects. 

At 7 (im, the situation is much different. Perhaps sur- 
prisingly, the PAH emission remains as good a star for- 
mation indicator in circumnuclear regions and starbursts 
as in spiral disks. The linear relationship previously de- 
fined by the disks alone holds up to S^Lyc ~ 10 4 L Q pc~ 2 ; 
adopting this value as transition point and assuming di- 



rect proportionality, we obtain: 

log(E 7/im ) = log(S Lyc ) - 0.25, < log(E Lyc ) < 4, (3) 

where again the size-normalized luminosities are in 
L pc~ 2 . The dispersion of circumnuclear regions and 
starbursts is 0.17 dex when galaxies beyond S^Lyc = 
10 4 L pc~ 2 are excluded, and increases to 0.25 dex (a 
factor 1.8) when the same objects as for the Sis^m^t - 
^Lyc relation are considered. The dispersion of disks is 
0.21 dex. Allowing the power-law index to vary, the least- 
squares fit over the entire < SLyc < 4 range gives an 
exponent of 0.96 ± 0.02. 

The 7 /im fluxes start to deviate significantly from 
the extrapolation of the linear correlation defined by 
disks above El yc ~ 10 4 L Q pc -2 . The starburst cores of 
NGC253, M82 and NGC6240 fall by a factor 2-3 be- 
low the expected values while Arp 220 lies more than an 
order of magnitude lower. Extinction effects alone can- 
not account for the saturation of the 7 (im emission be- 
yond SLyc = 10 4 Lq pc~ 2 , which is most probably caused 
by disappearance of the band carriers from the starburst 
cores (see Sect.[5j). This is not in contradiction with the 
different relation found at 15 /im: the fact that VSGs are 
larger and more resilient than PAHs allows the 15 /jm di- 
agnostic to continue rising up to higher star-formation ac- 
tivity levels than the 7 /an emission, though eventually 
VSG destruction might become significant too. 

Although extinction is not expected to be a dominant 
cause of the 7 /im emission deficit at high values of ^Lyc, it 
could, however, increase the dispersion at the highest star 
formation rate densities, together with variations in the 
average physical conditions of the gas (electronic densities 
and temperatures) and in metallicity. Another, certainly 
more important source of scatter is due to limitations of 
the available H line measurements and data used to es- 
timate the nebular extinction. In particular, although we 
tried to minimize such effects as much as possible, the 
apertures are not perfectly coincident, the angular resolu- 
tions are not perfectly matched, and the extinction correc- 
tion was sometimes derived in a region much smaller than 
the aperture. In view of all the uncertainties arising from 
use of inhomogeneous data and from assumptions about 
physical conditions in H II regions, and considering the fact 
that the observed dispersions are very small compared to 
the dynamic range of our relations, we have demonstrated 
that the two dust tracers investigated here constitute sat- 
isfactory and quantitative star formation estimators. We 
insist that the galaxies included in our sample are all of 
near-solar metallicities, while the relation between ioniz- 
ing photon luminosity and dust emission may be very sen- 
sitive to a decrease in carbon abundance. 

4.2. MIR color as an indicator of compactness 

Since the two dust tracers behave differently (the linear- 
ity ranges and dispersions of the relations discussed above 
are different), their variations relative to each other may 
provide useful diagnostics on the star formation activity. 
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Fig. 4. Evolution of the Li 5 ^ m ct / L 5 _ S 5 ^ ratio with in- 
creasing Lyman continuum luminosity per unit projected 
surface area. Apart from the disks, only galaxies with a 
CVF spectrum are shown here. 



Figure 0] shows how the Lis^m.ct/is-s.s^m ratio varies 
with increasing SLyc- It should be noted here that the un- 
certainty on the actual size of the emitting regions, which 
affects only the abscissa, is potentially quite large. It is 
however expected to be at most a factor of a few, i.e. 
very small compared with the variation amplitude of SLyc > 
which is three orders of magnitudes for the sole starburst 
regions. We find that the Sis^n^ct/Ss-s.s^m ratio, trac- 
ing to the first order the ratio of VSG emission to PAH 
emission, increases regularly from disks to mild starbursts 
to extreme st a rburs ts. A similar relationship was found 
by Dale et al in nuclear regions and extranuclear 

Hn regions of spiral galaxies. Their ionizing photon flux 
densities, derived in a uniform way from integral-field Pa/? 
and Br7 lines, correspond to ^Lyc, the quantity used here, 
between 10 2 ' 5 and 10 L Q pc~ 2 . Although with a large dis- 
persion, in part because of the inhomogcncous nature of 
the data used and the moderate angular resolution in the 
infrared, we show that this trend exists over a much larger 
range of star formation rate density. 

El Y c is essentially a quantification of the compact- 
ness of star formation activity, and is affected both by 
geometry effects (filling factor of the interstellar medium 
by Hn complexes) and by excitation effects (mass spec- 
trum of ionizing stars). These two effects produce qualita- 

ratio, the 



tively similar results on the £15 



5 /4m, ct/ £5— 8.5 |im 

PAH emission being reduced because the sites of emission 
(mainly the surface of molecular clouds) collectively de- 
crease in volume relative to H 11 regions, and because the 
PAH carriers may be destroyed over large spatial scales 
when the radiation field intensity becomes high enough. 
Excitation and d ensity effects are intimately connected 
( Dal e et al >ut variations in radiation hardness (or 

age of the dominant stellar populations) are also expected 



to incre ase the color scatter, as was found bv lRoussel et alJ 
l|2001bf) for central regions of spiral galaxies. The influ- 
ence of metallicity differences is expected to be negligible 
for this particular sample, and differential extinction ef- 
fects are very small compared with the dispersion, except 
perhaps for Arp220. 

5. Summary and conclusions 

At mid-infrared wavelengths, dust species emitting aro- 
matic bands seen mainly in the 6-13 /^m range and a con- 
tinuum rising toward longer wavelengths provide impor- 
tant observables for studies of star formation in dusty en- 
vironments, and provide finer details than the far-infrared 
emission of grains in thermal equilibrium, because of much 
more favorable angular resolution and lesser source confu- 
sion. The fraction of the total infrared power produced by 
each of these two species is of the order of 10-20 %, depend- 
ing o n the excitation condi tions of dust grains ( Dale et all 
l200ll iDale fc Helovl 120021) : the fraction contributed by 
very small grains, in particular, is remarkably constant 
up to very high average temperatures. In order for the 
SFR calibration at 5-8.5 /im derived here to be consistent 
with the calibration in terms of total infrared emission 
(between 8 and 1000 fim) proposed bv iKennicuttl l|l998h 
for starbursts, the power emitted in the 5-8.5 /^m range has 
to amount to 18% on average of the total infrared power 
(for galaxies less active than M82). Owing to the fact 
that the far-infrared emission of galaxies is not spatially 
resolved by IRAS, and that the central regions selected 
here emit only a fraction of the integrated MIR emission 
of each galaxy, we cannot rigorously estimate the part of 
aromatic bands in the energy budget separately for disks 
and for galactic centers (but this will become feasible in 
local galaxies with observations by the Spitzer satellite). 
We only remark that a power fractio n of 18% is somew hat 
larger than the fractions inferred bv lDale et alJ l|200lh for 
a wide range of dust temperatures. 

We have investigated the response of these dust species 
to the radiation field generated by massive stars, estimated 
independently and corrected for extinction, in a sample 
of star-forming sources of near-solar metallicity. In our 
sample, ionizing photon flux densities span a very wide 
range, from w 1 to « 10 5 5 L pc~ 2 . The regions consid- 
ered here are spiral disks on one hand, representative of 
quiescent environments, and circumnuclear regions on the 
other hand, extended on spatial scales of the order of the 
kiloparsec. 

Even though aromatic band carriers are on average 
heated by softer radiation than very small grains, we have 
shown that they can be used as a quantitative star for- 
mation tracer, their emission scaling linearly with the in- 
trinsic emission of hydrogen recombination lines over a 
dynamic range of four orders of magnitude in ionizing pho- 
ton flux densities. The relation found here confirms and 
extends that previously found for spiral disks up to much 
higher star formation rate densities. The global emission 
from aromatic bands starts to be damped past activity 
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levels only just milder than that of M 82. By analogy with 
what is observed in and around H II regions in the Galaxy 
and the Magellanic Clouds, this saturation is most prob- 
ably caused by the gradual destruction of aromatic band 
carriers effected by more and more intense far-u ltraviolet 
radiation fields <|TraJl998llContuTsi et aiTl2000tl In fact, 
this may be an indirect cause. iGiard et al.1 l|l994j) having 
found a tighter relationship of the 3.3 /im PAH brightness 
with the ionized gas density than with the radiation field 
intensity. Additional agents of dust grain destruction may 
be found in the enhanced c osmic ray density from numer- 
ous supernova explosions (Mcnncl la et al.l Il997f) . a nd in 
high- velocity starburst winds IjNormand et alJ ll995) . 

Such a behavior as seen here in galaxies was previ- 
ousl y reported for i ndivid ual photodissociation regions 
by iBoulanger et alJ |l998). The approximate threshold 
at which they observe a significant depletion of aromatic 
bands is ~ 10 3 5 times the radiation field of the solar neigh- 
borhood Go, or 10 4 ' 2 Lq pc -2 in t he unit used here (a dopt- 
ing Go = 2.2 x 1(T 6 WnT 2 from lMathis et alll983|) . The 
threshold applicable to galactic starburst regions, occur- 
ring arou nd IQ 4±0 - 5 L^pc~ 2 , is fully consistent with that 
found bv IBoulanger et all l|l998)) . It should be noted that 
resolution and dilution effects, as well as incomplete sam- 
pling of the explored radiation field range, hamper equally 
both studies, so that the actual value of the threshold 
is somewhat uncertain. The collective behavior of star- 
forming regions, integrated over kiloparsec scales, is nev- 
ertheless similar to that of individual H II regions and the 
associated neutral material surrounding them. This result 
suggests that the volume ratio of ionized regions on one 
hand, and surrounding regions where aromatic bands are 
excited on the other, does not vary in a systematic way 
up to the above mentioned radiation field intensity thresh- 
old, and then increases steadily, ionized regions occupying 
a growing fraction of the interstellar medium and starting 
to overlap. 

The continuum of very small grains (sampled at 
15 /im), on the other hand, provides a star formation rate 
tracer that is valid at higher radiation field intensities. 
In practice, variations of the spectral energy distribution 
of very small grains with their temperature distribution 
may cause appreciable deviations according to the sam- 
pled wavelength range, but we have shown here that the 
proportionality between ionizing photon fluxes and the 
15 fan continuum is impressively tight, as soon as the 
VSG continuum dominates the bandpass, and at least up 
to ELyc = 10 5 L©pc -2 . Very small grai ns may also be 
destr oyed in very harsh radiation fields ijContursi et alJ 
2000), but this effect is not observed here except possibly 
in Arp220, where it is however not separable from optical 
depth effects. 

New space missions such as Spitzer are making the 
mid-infrared window ever more accessible and are going 
to perform large surveys of galaxies. The choice to measure 
the continuum of very small grains at 15 (im was dictated 
by the limited wavelength coverage of the data we used. 
However, with the Spitzer satellite, this continuum will 



be observable primarily through the 20-28 /im filter of the 
MIPS instrument, and through the IRS Long-Low spec- 
trometer for brighter galaxies. For local galaxies, these 
wavelengths promise an excellent star formation tracer, 
following a single regime from disk-like to very high radi- 
ation field intensities. The MIPS 24 /im filter will detect 
the continuum of very small grains from z = to z « 0.6, 
shifting gradually down to 15 /im, then the aromatic band 
cluster at 6-9 /im at z w 1.8-2.7. The MIPS 70 fan filter 
will cover the continuum of very small grains from 30 /im 
to 15 /im from z w 1.3 to z w 3.5. The quantitative re- 
lationships that we have derived in this paper might thus 
prove very useful in the immediate future. 
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Appendix A: Notes on individual sources and 
photometry 

We give here additional details concerning individual 
sources as well as MIR and line flux measurements, for 
the spiral galaxies whose circumnuclear regions are stud- 
ied in this work and the starburst systems. Whenever 
the Ha data included a contribution by the adjacent 
[Nil] AA 6548, 6583 A lines, we applied a correction factor 
of 0.75 to obtain pure Ha fluxes (see Sect. E2J>. "[Nil]" 
refers to both lines at A = 6548 and 6583 A, except if 
the wavelength of the line actually meant is given. Unless 
explicitely specified otherwise, extinction values derived 
from H line measurements are referred to a uniform fore- 
ground screen (UFS) model, case B recombination, with 
n c = 100 cm" 3 , and T c = 5000 K. 

We brought the MIR maps in all bandpasses to 
the same angular resolution before performing aperture 
photometry, substituting the extended non-gaussian PSF 
with a gaussian PSF of FWHM 6" (when the pixel size 
is 6"), 3.5" (when the pixel size is 3") or 3" (when 
the pixel size is 1.5"). To do this, we used an iterative 
procedure with a gain of 5% to ensure convergence, 
centered on the brightest pixel at each step. When using 
a map to measure hydrogen recombination line fluxes, we 
then convolved this map to the same angular resolution 
as in the MIR. Except when contrary indication is 
given below, we chose an homogeneous definition of the 
aperture as 2.5 times the half-power beam width fitted 
on the central MIR brightness profiles; this provides 
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a very good approximation of the total flux of the 
central regions, coin ciding well with sizes reported in 
iRoussel et all l|2001rJ) (obtained by decomposing bright- 
ness profiles into a gaussian core and an exponential disk). 

NGC986 - We assumed that the extinction derived 
from the H/3/Ha ra t io in the central 4" x 4" from 
IVeron-Cettv fc VerorJ (|l986h is re presentative of that in 
our la rger aperture. The data of IVeron-Cettv fc VerorJ 
ljl98fih imply [Nil] A 6583 A/ (Ha + [Nil] A 6583 A) =0.30 
at the nucleus. 

NGC 1097 - This strongly barred Liner/Seyfert galaxy 
has a bright star-forming ring of diamete r ft 20" (e.g. 
iHummel et al.llT987tlKotilainen et al.ll200o(l . The nucleus, 
which is resolved and separable from the ring in our maps, 
cont ributes negligibly to the integrated Ha and Br 7 emis - 
sion l|Storchi-Bergman et al.ll996tlKotilainen et al.l2000j) . 
as well as to the total MIR emission. We used an aperture 
of 45", encompassing the whole emission from the ring 
l|Roussel et al.ll2f)01bl) . We corre cted for the aver a ge ex- 
tinction based on the results of iKotilainen et aJl J2000) 
derived from Ha/Br 7 ratios. We recomputed the weight- 
ing by the Br7 luminosity, with adjustments for the [N 11] 
contribution to their Ha data and the different extinction 
laws adopted. 

NGC 1365 — We combined the Br7 measurement of 
IPuxlev et al.1 (|l98St) with the Ha flux integrated within 
the same region from an Ha + [N 11] map to derive the 
extinction, and assumed that it represents accurately the 
extinction within our larger aperture. The Seyfert nucleus 
does not contribute importantly to the H line and MIR 
dust emission. The total Ha flux in the central 4" x 4" 
ijVeron-Cettv fc VerorJ Il98^) . which also includes emis- 
sion from adjacent "hot spots," is only 8% of that in 40". 
MIR diagnostic line ratios suggest that star formation ac- 
tivity dominates the low excitation (< 50 eV) line spec- 
trum at these wavelengths as w ell as the MIR an d far- 
infrared continuum luminosities 1 Sturm et alJl200^) . The 
nucleus is unresolved in the ISOCAM maps, preventing 
an accurate estimate of its contribution to the MIR fluxes, 
but the ISOCAM CVF data do not provide evidence for 
a significant AGN contrib utio n based on the diagn ostics 
of lRigopoulou et alJ <|l999^ and lLaurent et all l)200(J) . The 
[Nil] /(Ha + [Nil]) ratio in the central 4" x 4" and in sev- 
eral hot spots within 1 4" x 20" is ft 0.3 jAlloin et alll98ll 
IVeron-Cettv fc VeroullTgsi . 

NGC 4102 - The nucleus is saturated in both ISOCAM 
LW2 and LW3 observations, more severely for LW2. This 
galaxy generat es a powerful outflo w detected in the Pa/3 
and Br 7 lines l|Roussel et al.ll2003l) . 

NGC '4293 - Since the central MIR source is very 
small compared to th e pixel size of the ISOCAM maps 
ijRoussel et al.lE001a|) . we do not use the fitted aperture 
for area normalization, but instea d the size derived from 
the Pa a and Ha maps. The data of lVeron-Cettv fc VerorJ 
in the central 4" x 4" give a high line ratio 
[Nil] A6583A/(Ha + [Nil] A 6583 A) = 0.68 due to the 



Liner nucleus, and an Ha flux accounting for ft 9% of the 
flux in d = 32.4". 

NGC 4691 - We combined the Br 7 measurement of 
IPuxlev et al. I (ll99Cfl with the Ha flux integrated within 
the same region from an Ha + [N 11] map to derive the ex- 
tinction, and assumed that it remains the same within our 
larger aperture. The [N i i] / (Ha + [ N III) ratio in the cen- 
tral hot spo ts is ft 0.3 l|Keellll983r iGarcia-Barreto et alJ 
Il995l Il999|) . As the central structure contains multiple 
knots which are partly blended in the ISOC AM maps, and 
is not well represented by a single gaussian l|Roussel et alJ 
2001a), we strayed from our general definition to deter- 
mine the area normalization. As the MIR and Ha emission 
of NGC 4691 lacks in the disk and is very diffuse outside 
the central star-forming knots, we simply selected pixels 
above the 3<r brightness level in the Ha map and added 
their areas to compute an equivalent diameter. 

NGC 5194 - We considered the central emission 
plateau of diameter 90". The motivation for this choice in- 
stead of selecting a smaller region around the nucleus was 
that the plateau represents a transition between disks and 
more active circumnuclear regions in terms of MIR prop- 
erties. Besides, the emission from the Seyfert nucleus is 
diluted and completely negligible within this large aper- 
ture. We adopted the average extin ction, weighted by in- 
trinsic Ha luminosities, derived by IScoville et alJ l|200lfl 
from Ha/Paa decrements of a large sample of Hn re- 
gions, assuming it is representative of the effective extinc- 
tion throughout the central plateau. We applied a small 
correction to this extinction to account for the different 
line emissivities and extinction laws adopted. The [Nil] 
/ (Ha + [Nil]) ratio in the central d <, 20" is high and 
reaches ft 0.85 due to the Sey fert nucleus but goes down 
to ft 0.3 outside these regions llRose fc Searlelll982ft . 

NGC 5236 - lOenzel et all lll998l) report in the SWS 
14" x 20" aperture an extinction of 5 mag, larger than 
the value adopted here. The ratio [Nil] A 658 3 A /(Ha + 
[Nil] A 6583 A) ft 0.3 in th e central ft 5" ijKeell Il984l 
IVeron-Cettv fc Veronlll98fJ) . The nucleus is saturated in 
the ISOCAM LW2 and LW3 maps, more severely for LW3. 
We thus used maps simulated from the CVF spectral cube 
to measure the LW2 and LW3 fluxes. We estimate in this 
way that the missing flux fractions due to saturation are 
11% and 34% respectively, assuming that the effects dis- 
cussed in Sect. l3.il making photometry from CVFs uncer- 
tain by 10-20%, are negligible. 

NG C 6946 - We combined the Br 7 flux of lPuxlev et alJ 
l|l988|) with the Ha flux integrated in th e same aper- 
ture from an Ha + [Nil] map. iKeell l)l984j) gives 
[Nil] A6583A/(Ha + [Nil] A 6583 A) ft 0.37 in the central 
8.1". We assumed the same ratio throughout our larger 
aperture to calibrate the Ha + [Nil] map. The nucleus 
is saturated in the ISOCAM LW2 and LW3 maps, more 
severely for LW2. As for NGC 5236, we used maps simu- 
lated from the CVF spectral cube, and we estimate that 
the missing flux fractions due to saturation are 27% and 
13%, respectively. 
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NGC7552 - IVerma et ail ll2003h published Bra and 
Br/? fluxes obtained with SWS in a 14" x 20" aperture 
which matches fairly well our circular aperture of 21.8". 
We used Bra only because of possible blending of Br/3 with 
H 2 1 — 0(2), and combined it with the Ha flux integrated 
over the same region from an Ha + [N n] map to derive the 
extinction. The nucleus is slightly saturated in the L W3 
map. From the data of IVeron-Cettv fc Veronl <jl986() for 
the central 4" x 4", [Nil] A 6583 A/ (Ha + [Nil] A 6583 A) 
= 0.37 and the Ha flux is 24% of the flux in 21.8". 

NGC7771 - Br7 and radio imaging shows that the 
central d w 10" area hosts the active star-forming re- 
gions, distributed mainly along a circumn uclear ring 
llNeff fc HutchingJl992HReunanen et alJ2000(l . Note that 
although iReunanen et al .1 (1200011 mention that they cor- 



rected for velocity shifting of the Br7 line outside the 
narrow -band filter they used, the data from iDale et all 
(2003), which we adopted, yield a Br7 flux almost three 
times higher in the central 5" x 10". 

NGC253 - The Sevfe rt nature of the nucleus 
I Veron-Cettv fc Veronll2001^ is unconfirmed by other op- 
tical spectroscopic studies, and by near- and mid-infrared 
spect roscopy (e.g. lEngelbracht et all Il998t ISturm et alJ 
2000j); neith er by our ISOCAM dat a bas ed on the di- 
agnos tics of iRlgQDOulou et all 1 19991) and ILaurent et alJ 
<|2000h . We used the Br 7 flux of lEngelbracht et all l|l998) 
integrated over d = 15", which contains all the flux and 
coincides with the fitted size of the MIR emission. We 
used their Pa/3/Br7 flux ratio measured in 2.4" x 12" 
to derive the extinction. Published extinction and Qhyc 
estimates from H line data vary greatly. For instance, 
IVerma et alJ l)2003l) derived from SWS line observations 
Ay IX ^ g maKi considerably lo wer than IX = 30 mag 
reported bv lGenzel et al.l ljl998|) . although the QLyc values 
are similar. 

NGC520 - Br7 line emission in this interacting sys- 
tem is detect ed at the primary nucleus within a r egion 
of ~ 5" x 3" (IStanfordlll99lUKotilainen et alJl200ll) . We 
combined the Br7 flux of IStanfordTl|l99lj) in the central 
6" x 8" with the Ha flux measured within the same aper- 
ture from an Ha+[Nn] map, to derive the extinction. We 
made use of the Hq /(Ha+[Nii]) ratio of 0.58 found by 
IVeilleux et all l)l995|) at the nucleus. Since the primary 
nucleus suffers from very large extinction, the use of Ha 
photometry in a larger aperture is especiall y uncertain. 

NG C1808 - We used the Br 7 flux of iKra.hhe etlil 
l(l994h integrated in a 20" aperture centered at the 
nucleus, which contains the quasi-totality of the flux 
and coincides reason ably well with the fitted MIR size. 
iKrabbe et alJ l)l994[) derived extinction values ranging 
from 3 to 5 mag from Ha/Br7 decrements, in excellent 
agreement with estimates from H/3/Ha decrements quoted 
in the same paper, so we applied an average of 4 mag. 

NGC3034 (M 82) - We used directly t he ex tinction 
and QLyc results of lForster Schreiber et al.l ( 200 ll) for the 
starburst core within d = 30". These results were de- 
rived from an extensive set of H lines from optical to radio 
wavelengths which are best fitted by a mixed model and 



with deviations from the iDraind l|l989|) extinction l aw, at 
A = 3 — 10 /Ltm, as found towards the Galactic Center ijLutd 
19993). The MIR source is elongated along the optical ma- 
jor axis. Computing the quadratic mean of the major axis 
and minor axis widths, we find that 2.5 x HPBW = 29.95" 
at 15 /xm, thus extremely close to the adopted aperture 
of 30". The 5-8.5 /im emission is more extended, and we 
would have derived an aperture of 35" from it. 

IC342 - We adopted an aperture of 17" x 17" cor- 
responding t o the size of the Br7 line map provided by 
Bokcr ct^y l)l997|) , which encompasses the circumnuclear 
starburst ring. This aperture is close to the intrinsic size of 
the mid-IR source fitted on the surface brightness profiles, 
between 19" and 22". The areas agree to within 30%, and 
the MIR fluxes to within 9% in the different apertures. We 
combined t he Br^ flux integra ted over the map with the 
Bra flux of IVerma etafl J2003) obtained in the 14" x 20" 
SWS beam (excluding their Br/3 measurement because of 
possible contribution from the H 2 1 — 0(2) line and Pfa 
because of larger uncertainties on the extinction law near 
7 fim) . Fits assuming a UFS and a mixed model are both 
well constrained, and the derived extinctions imply nearly 
identical QLyc (within 1%). 

NGC3256 - The MIR emission can be separated into 
two components of different sizes, a core of w 8" super- 
posed onto a source of « 20". Therefore, a single gaus- 
sian does not provide a good fit to the surface bright- 
ness profile. We adopt the size o f the larger source as 
our aperture. iKawara et al ] lll987t) report a Br7 flux of 
14.8 x 10 -17 Wm" 2 in a smaller aperture of 9" x 18", al- 
most identic al to the total BrT flux of 15 x 10~ 17 Wm~ 2 
measured bv lMoorwood fc Olival (^994) in a field of view 
of 34" x 34". We derive the extinc tion in the c entra l 
3.5" x 3.5" from the Pa/3/Br7 ratio of iDovon et aD l)l994h 
and assume that it does not vary inside our larger aper- 
ture. 

NGC6240 - The MIR emission of NGC 6240 is unre- 
solved in the ISOC AM data. W e adop t a size of 3" based 
on the Br7 map of iTecza et all ll2000h . As a comparison, 
we derive from brightness profile fitting a HPBW of the 
order of 2.7", which is ill-constrained since it is smaller 
than t he cleaned PSF HP BW of 3". We adopted the Br 7 
flux of iRieke et alJ lll985l) measured in a 8.7" aperture, 
and the Pa/3 flux of ISimpson et alJ l)l996j) measured in 
a slit of width 1.5" oriented al ong the major axis of the 
source. In view of the Br7 map of lTecza et al.ll(2000l) . these 
apertures should in clude the total fluxes. The Paa flux of 
Ide Pov et alJ l|l986h . measured in a 5.5" aperture, yields 
a negative extinction when combined with the Pa/3 flux; 
it may thus be strongly underestimated. MIR diagnostics 
indicate that starburst activity dominates the dust emis- 
sion and the low-excitation fine-structure line emission at 



these wavelengths ( 


Genzel et a,l.||l998l iRlgopoulon et alJ 


1999: ILaurent et al. 


2000;: but see alsolLutz et al.ll2003). 



IRAS 23128-5919 - The two galaxies in this merging 
pair are separated by a projected distance of less than 
5", and are marginally resolved by ISOCAM. We derive 
HPBWs of the order of 2.5" for both nuclei. For lack 
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of adequate high-resolution data, and in analogy with 
NGC 6240, we adopt a total size of 3" for each nucleus, 
thus an equivalent size of 4.2". To derive the extinction, 
we combined measurements of pure Ha and Paa fluxes 
in the southern nucleus, which is much brighter than the 
norther n nucleus. Note that because of the narrow slit 
used bv lDuc et alJ l|l997l) . the Ha fluxes may be underes- 
timated. In the MIR, we performed the brightness profile 
fitting and aperture photometry of the two blended nu- 
clei together. Since the CVF spectrum has poor signal to 
noise ratio, and the 12-18 (im flux simulated from the CVF 
is not in agreement with the observed broadband flux, 
we used only broadband images. From both optical and 
MIR diagnostics, this system is classified as a starburst 
in which the large velocities observed in emission lines of 
the southern nucleus are caused by starburst s uperwinds 
(Johansson & Bergvalflll988l:lLutz et al.lll999bl) . 

Arp 220 - The two nuclei of this merger are separated 
by ~ 1" and are totally unresolved in the ISOCAM data. 
The size that we derive from 15 fim brightness profile fit- 
ting is 1.8-1.9" HPBW. We use a total extent of 2", as 
implied by the high angular resolution observations of the 
mid-IR emission bv lSoifer et al.1 l)l999|) . MIR diagnostics 
indicate that starburst activity dominates the dust emis- 
sion and the low-excitation fine-structure line emission 
llGerrzel et al!l998tlRigODOulou et alll999ULaurerit et all 
|200(A see also the detailed study bv lSpoon et alll2003|) . 
The extinction that we derive from the Bra and Br7 
fluxes is in good agreemen t with that predicted by 
Anantharamaiah et al.l l|2000l) from a set of radio recom- 
bination lines. 



Appendix B: Spectral observations of IC342 

The CVF observations of IC 342 consist of a single point- 
ing centered on the nucleus, with the 3" pixel -1 scale giv- 
ing a 96" x 96" field of view, and contain w 35 exposures of 
2.1s per wavelength channel for a total integration time of 
3 h 18 m . A representative spectrum, in a region of 21" x 21", 
is shown in Fig. IB. II It shows all the usual features of star- 
forming galaxies, with aromatic bands nearly invariant in 
shape and relative intensity (see Fig. Notice the faint 
emission bands that are also present in the spectra of the 
other galaxies of this sample, and which are attributable 
to carriers of the same nature as those emitting the bright 
bands at 6.2, 7.7, 8.6, 11.3 and 12.7/xm. Those contribut- 
ing to the emission within the LW3 filter and observed in 
this spectrum are d etected at 12.0, 13.6, 14.3 and 15.7/xm 
l|Sturm et al.ll2000j) . 
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